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Abstract 


This paper highlights the Advanced Rocket Injec- 
tor /Combustor Code (ARICC) that has been developed at 
Rocketdyne to model the complete chemical/fluid/ 
thermal processes occurring inside rocket combustion 
chambers. The code, derived from the CONCHAS-SPRAY 
code originally developed at Los Alamos National Labor- 
atory (Ref. 1) incorporates powerful features such as 
the ability to model complex injector combustion cham- 
ber geometries, Lagrangian tracking of droplets, full 
chemical equilibrium and kinetic reactions for multi- 
ple species, a fractional volume of fluid (VOF) 
description of liquid jet injection in addition to the 
gaseous phase fluid dynamics, and turbulent mass, 
energy, and momentum transport. Atomization and drop- 
let dynamic models from earlier generation codes are 
transplanted into the present code. Currently, ARICC 
is specialized for liquid oxygen/hydrogen propellants, 
although other fuel/oxidizer pairs can be easily sub- 
stituted. ARICC is a two-dimensional/axisymmetric 
real-time transient code, and effort on developing a 
three-dimensional counterpart is underway. 

Introduction 


Combustion processes occurring in a fluid medium 
are probably among the most complex of all flow phen- 
omena. The complexity stems from two primary reasons. 
Firstly, the so-called "combustion" process actually 
encompasses many physical processes of different types. 
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different natures, different temporal and spatial 
scales, and different degrees of being describable by 
deterministic models. Secondly, all of these processes 
are usually strongly coupled together, making it dif- 
ficult to reduce the problem down to its individual 
component processes. Any attempt to numerically model 
such complex combustion phenomena would need to make 
careful distinction between those processes that can 
be considered as the primary driving ones, and hence, 
should be modeled as vigorously as possible, and those 
that can be considered secondary and are included only 
for their usefulness in "interpreting" the situation 
into the desired results. Such decisions are often as 
important as the quality of the mathematical models 
themselves . 

Inside a liquid rocket engine, the process usually 
starts with the injection of fuel and oxidizer into the 
combustion chamber in both liquid and gaseous states 
under high pressure and at greatly different velocities 
so as to encourage rapid atomization of the liquid com- 
ponent. The flow created by such injection is highly 
turbulent and rotational. Furthermore, in the case of 
liquid oxygen/hydrogen engines, the propellants are at 
temperature extremes ranging from cryogenic (70 K) to 
superhot (3000+ K) . Thus, even in the precombustion 
state the flow is fully elliptic and composed of multi- 
ple species with real gas /liquid properties. The in- 
ception of combustion magnifies the local temperature, 
pressure, and species concentration gradients, and adds 
one more time scale to those of the bulk flow tran- 
science and of the turbulence. In light of these cir- 
cumstances, a list of requirements deemed necessary for 
the numerical modeling of the combustion process has 
been generated and are listed in Table 1. 

The ARICC code represents an attempt to build a 
modular model to describe the complete combustion 
process inside a liquid rocket engine using state-of- 
the-art finite difference techniques. It is an ad- 
vancement over earlier-generation performance codes, 
many of which are used as industrial standards, that 
frequently describe one aspect of the combustion 
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Table 1. Code Features Required for 
Liquid Rocket Engine Modeling 


NECESSARY FEATURES 

FEATURES INCORPORATED INTO ARICC 

• VISCOUS ELLIPTIC FLOW 

• FULL NAVIER-STOKES EQN. + SGS TURBULENCE 
MODEL, ALL FLOW SPEEDS. 

© FULL CHEMISTRY 

© MULTIPLE SPECIES KINETIC AND 
EQUILIBRIUM REACTIONS 

® LIQUID JET INJECTION/ 
2 PHASE FLOW 

® VOF TWO-PHASE, FREE SURFACE 
DESCRIPTION 

« DROPLET DYNAMICS 

« MONTE CARLO SCHEME/LAGRANGIAN DROPLET 
TRACKING. DROPLET HEATUP 
AND VAPORIZATION 

© ATOMIZATION 

» CICM ATOMIZATION MODEL 

• FINE RESOLUTION IN 
CRITICAL AREAS 

© ARBITRARY GRID WITH CLUSTERING. 


process in great detail (e.g., chemistry) but make 
gross assumptions about the other aspects (e.g., 
neglecting flow recirculation by using the stream tube 
method). These codes (Ref. 2-6) with their primary 
limitations are listed in Table 2. However, they all 
contain valuable experience and empirical models that 
have and will continue to be assimilated into the 
ARICC code. 


Table 2. Earlier-Generation Rocket 
Combustion Codes 


| CODE NAME 

PRIMARY FUNCTIONS 

MAJOR LIMITATIONS 

m 

STEADY-STATE ROCKET COMBUSTION 
OF H 2 /LOX FOR COAXIAL INJECTION 

STREAM TUBE FLUID DYNAMICS. 
ATOMIZATION CORRELATIONS FROM 
J-2 ENGINE 

TRANS78 3 

MULTIPURPOSE CHEMISTRY CODE FOR 
ROCKET PERFORMANCE, DETONATION, 
CONSTANT P OR CONSTANT V 
COMBUSTION ETC. 

ID EQUILIBRIUM ISENTROPIC. 

NO SPATIAL RESOLUTION OR FLOW 
DYNAMICS 

SDER 4 

STEADY-STATE ID FLUID DYNAMICS 
WITH 3D DROPLET DYNAMICS FOR 
IMPINGING TYPE ELEMENTS. LIQUID/ 
GAS, SUB. AND SUPER-CRITICAL 
VAPORIZATION OF DROPLETS 

DROPLET CORRELATIONS DEVELOPED 
FROM COLD FLOW. STREAM TUBE 
FLUID DYNAMICS 

T pp5 

TIME-DEPENDENT PROGRAM FOR 
SMALL PULSED TYPE ENGINES WITH 
HYPERGOLIC PROPELLANTS. 

3D DROPLET DYNAMICS. 

ID FLUID DYNAMICS. IMPINGING 
LIQUID-LIQUID ELEMENTS ONLY. 
SUBCRITICAL DROPLETS 

GKAPS 

STEADY-STATE CHEMICALLY REACTING 
GAS MIXTURES OF UP TO 150 SPECIES 
AND 70 NONEQUILIBRIUM REACTIONS. 

ID FLOW. 

ADIABATIC COMBUSTION CHAMBER. 
NO DROPLETS. 
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Process of Code Development 


The process of code development consists of four 
steps: (1) identify the critical processes; (2) de- 

termine which of these processes can be described 
vigorously and which should be modeled empirically, 
and devise mathematical expressions for each category; 
(3) translate these expressions into numerical formu- 
lations and assemble a logical algorithm for evaluating 
them and displaying the results; and (4) design steps 
for verification of the model components and for 
"anchoring" or fine-tuning of the empirical model 
constants . 

The above process is depicted in Fig. 1. As a re- 
sult of the initial scrutiny in step one, several 
assertions can be made concerning liquid rocket combus- 
tion. First, of the two primary transport processes of 
diffusion and convection (radiation is not important at 
these temperatures), the former controls the transfer 
of species with respect to each other while the latter 
controls the bulk transfer of mass, momentum, and 
energy. As the liquid oxygen/hydrogen reactions are 
relatively fast, the steady-state flame is thus diffu- 
sion controlled and can be described sufficiently with 
equilibrium reactions. Viscosity (diffusive transfer 


OBSERVATION INTERPRETATION CONCEPTUALIZATION DISCRETIZATION 



Fig. 1. Steps in Formulation of a 
Computer Model 
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of momentum) is important only in stagnant zones of the 
combustion chamber where recirculation controls wall 
heat transfer, and indirectly through the turbulence 
scale which affects the species diffusion and atomiza- 
tion processes. Second, the overall performance of the 
rocket engine is strongly dependent on the effective- 
ness of the atomization process . Unvaporized and hence 
uncombusted droplets are the greatest cause of loss of 
engine efficiency. Third, while the steady-state flame 
is adequately described by a diffusion-controlled equil- 
ibrium flame model, the ignition and flame holding mech- 
anisms and wall-streaking phenomenon of concern can only 
be understood with detailed consideration to the local 
geometry and flow patterns and transient heat transfer 
processes. The life and durability of the engine de- 
pends largely on these last mechanisms. 

In other words, an adequate numerical model must 
have the capability to fully describe the following 
types of phenomena: viscous flow as described by the 

full Navier-Stokes equations with turbulence modeling; 
multiple species chemical kinetics, liquid jet injec- 
tion, atomization, and subsequent vaporization and 
transport of the droplets; and fine grid resolution 
near boundaries and corners of particular interest. 

All of these capabilities have been built into ARICC. 

The type of mathematical model used for each phenomena 
is listed in Table 3. 

Before a more detailed description is given of the 
code, a comment about the code's verification is in 
order. The basic fluid dynamics, which forms the basis 
of the solution scheme, will and is being checked out 
with experimental laser velocimetry measurements done 
on a low pressure cold flow test chamber under the cur- 
rent contract, and the final results will be reported 
in a future writeup. Optical Schlieren and infrared 
data obtained in the same tests also provide a qualita- 
tive means of verifying the overall temperature pro- 
files. Far more critical however are the verification 
needs pertaining to the empirical models , where the 
greatest uncertainty lies and where the greatest varia- 
tion from problem to problem still exists in the form 
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Table 3. Mathematical Formulation Used for 
Various Physical Processes 


1. FLUID THERMODYNAMICS 


MASS CONSERVATION (CONTINUITY) 

PD 

MOMENTUM CONSERVATION (NAVIER-STOKES) 

PD 

ENERGY CONSERVATION 

PD 

II. PHYSICAL MODELS 


EQN. OF STATE (CONSTITUTIVE) 

AL/TA 

TURBULENCE 

AL/PD 

CHEMISTRY 


RATE EQUATIONS 

AL 

EQUILIBRIUM EQUATIONS 

AL/TA 

RADIATION 

PD 

III. DROPLET DYNAMICS 

OD 

DISCRETE FORMULATION 

OD 

CONTINUUM FORMULATION 

PD 

TYPE OF GOVERNING EQN. 


PD = PARTIAL DIFFERENTIAL 


AL = ALGEBRAIC 


TA = TABULAR 


OD = ORDINARY DIFFERENTIAL 



of arbitrary coefficients. Replacing these coeffi- 
cients with universal constants or expressions is a 
task of utmost concern. This applies to the turbu- 
lence model, supercritical droplet properties, the 
atomization model, as well as the droplet heatup and 
vaporization models. Few correlations are available 
and measurements are difficult to make and. hard to 
come by, especially in the area of droplet dynamics. 
Finally, the choice of chemical equations, species, 
and rate or equilibrium constants to describe the 
chemical processes also represents a major area of 
uncertainty. In the case of oxygen/hydrogen kinetics, 
the role of the third body reactions is largely un- 
determined. An initial one-dimensional test case 
using the ARICC code however does indicate that the 
resulting temperature profile inside a premixed com- 
busting tube is strongly influenced by the omission 
of certain third-body reactions. Much study is still 
required before any conclusion can be made. 
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Code Description 


Code Organization 

Figure 2 shows the organization of ARICC, which is 
typical of a full-fledged combustion code. As pre- 
viously stated solution of the fluid dynamic governing 
equations of mass, momentum, and energy with attending 
boundary conditions form the core of the numerical 
algorithm. On top of that, the empirical models 
describing the gas or liquid properties (enthalpy 
tables, vapor pressure curve fits, etc.), the chemical 
reactions, the atomization and droplet heatup and evap- 
orization processes, and the turbulence model are called 



Fig. 2. Current ARICC Code Organization 


as subroutines as and when they are required. To make 
an operationally useful code, however, the code must 
be given a flexible input /output capability for grid 
generation, field initialization, and graphical and 
tabular display of outputs. All these functions are 
in turn controlled by a master main program which also 
takes the responsibility of ensuring stability through 
time-step size control and various artificial numeri- 
cal smoothing schemes. Figure 3 is a simplified flow 
chart of ARICC. 
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Fig. 3. Simplified Flow Chart ARICC Code 


Basic Solution Algorithm 

The basic Implicit Continuous Eulerian with Arbi- 
trary Lagrangian-Eulerian (ICE’d-ALE) mesh treatment 
solution technique used in ARICC has been described in 
detail in Ref. 1 and will not be repeated here. The 
major complication involved in the ARICC code is the 
need to handle both compressible gases and incompres- 
sible liquid within the numerical framework. To do 
this, cells that contain liquid must first of all be 
distinguished from cells that do not. To this end, 
the Fractional Volume of Fluid method of description 
has been adopted from the free surface flow SOLA-VOF 
program (Ref. 7). A cell flag variable (F) is used to 
represent the fraction of the volume of a cell that is 
occupied by compressible gases. Thus, three differ- 
ent kinds of active cells are possible (Fig. 4). 

Cells that are totally or partially gaseous (with F 
values greater than 0 and less than or equal to 1) 
will be treated as a "normal" cell using the standard 
ICE'd pressure iteration scheme. However, the cell 
volume used in the calculation will be that of the 
effective gaseous volume equal to F times the total 
active cell volume. Cells with flag value equal to 
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• ACTIVE CELLS DISTINGUISHED BY THEIR FLAG VALUES: 



FULLY LIQUID 
F = 0.0 


PARTIALLY LIQUID 
F = 0.7 


FULLY GAS 
F = 1.0 


DEF F = 


NET GASEOUS VOLUME (EXCL. DROPLETS) 
LIQUID VOL + NET GAS VOL. 


• DROPLETS MAY EXIST IN CELLS WITH F > 0.0 

• SINGLE-VALUED THERMODYNAMIC PROPERTIES FOR ALL CELLS 

• PRESSURE IN FULLY LIQUID CELLS DETERMINED BY PRESSURE 
ITERATION SCHEME FOR IMCOMPRESSIBLE FLUID 


Fig. 4. Different Types of Active Cells 
in Two-Phase ARICC 


0 are fully liquid cells, and while they are treated 
the same way in the rest of the program, in the pres- 
sure iteration step these cells are given a "pseudo- 
incompressible" treatment by using a pressure iteration 
scheme that tries to enforce a conservation of the 
fluid volume, and hence density (Ref. 8). No equation 
of state relation is applied to the liquid portions and 
a constant density (as well as temperature) is simply 
assumed. Chemical reactions are permitted only inside 
the gaseous phase and no direct mass or heat transfer 
takes place across the gas /liquid interface. Thus, 
other than for the fluid dynamical interactions, the 
only path by which the liquid affects the gaseous mix- 
ture is by first turning into droplets and evaporating. 
The atomization model is described below. 

Since the main reason for inclusion of the liquid 
phase description is for modeling of the liquid jet 
column coming in through the coaxial injector of a 
rocket engine, the bulk of the liquid is currently 
assumed to be always around the central axis of the 
model. This simplifies the determination of the orien- 
tation of the liquid-gas interface necessary for the 
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calculation of the flux of liquid from one partially 
liquid cell to the next, and also simplifies the 
graphics scheme for plotting the interface. A gener- 
alization to accommodate the most arbitrary liquid- 
gas combination can be readily made. 

Atomization Model 

At designated intervals, numerical droplet groups 
representing a number of particles of the same proper- 
ties are created from the cells that contain a liquid- 
gas boundary. The atomization model consists of cal- 
culation of the stripping rate and the mean drop size 
(and possibly a statistical distribution function of 
the drop sizes) from the local fluid dynamic proper- 
ties. Currently, the model used in ARICC is that 
taken from the CICM program (Ref. 2). The expressions 
are as follows: 


Mean 

Stripping Rate 


m a = c a 


2 2 

Wr > 

y p j 


ll/3 


ttD^. (AZ) 


Mean 

Drop Size 


D. = B a 

3 A 


pu 2 

g r 


1/2 


2/3 


where 


c a’ b a 


u 


V P 3 

P 3 


empirical constants 
jet diameter 

gas velocity relative to liquid jet 

axial location 

gas, jet density 

jet viscosity 

jet surface tension 


At creation, the droplets are given a temperature 
and density equal to those of the liquid from which 
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they are created, and an appropriate amount of liquid 
is eliminated from the cell liquid volume. The re- 
verse process of droplet "recondensation" into the 
liquid is also allowable, if the droplet velocity 
should lead it to head back toward a liquid surface. 
Also, the droplets are given an initial speed equal to 
that of the local liquid surface and a direction that 
randomly spans a 90-degree fan, ranging from being 
aligned with the liquid surface to being radially 
away. Obviously, these fine details of the model are 
somewhat arbitrary and can be changed as more exper- 
imental data becomes available. 

Both because of the fact that the complete atom- 
ization process from liquid medium is being modeled 
and because of the substantial density and sizes of 
droplets dealt with in liquid rocket engines (up to 
100 microns as compared with submicron droplets 
usually encountered in gas turbines) , the droplets in 
ARICC are given a finite volume that would physically 
exclude gas and liquid from within the cell. Dynam- 
ically speaking, the "bouyancy" term is being in- 
cluded. Thus, the droplet volumes enter into the de- 
termination of the effective (gaseous) cell volume. 
This situation is also depicted in Fig. 4. However, 
a droplet can only be in one computational cell at 
any time, as its location is uniquely defined by a 
single set of coordinates. Currently, interaction 
among droplets (e.g., coagulation) and film forma- 
tion along solid surfaces as the droplets hit the 
walls are not modeled. 

This advanced description of the droplet spray 
plus the two-phase liquid jet description and the 
CICM atomization model make the ARICC code unique 
in its capability for modeling liquid rocket engines. 
The entire combustion process can now be analyzed, 
in whole rather than in parts. It combines the dif- 
ferent branches of fluid dynamics, namely, two-phase 
flow, reactive flow, particulate flow, and compres- 
sible recirculating flow into a single transient 
finite difference model. The atomization process 
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that is so important to rocket engine performance can 
now be studied. Finally, the droplet finite volume 
displacement and recondensation features pave the way 
for simulation of droplet breakup and coagulation 
phenomena. 

C::ygen/Hydrogen Chemistry Model 

Although the primary reactions of the oxygen/ 
hydrogen chemistry are well known, the secondary 
reactions, especially those involving third bodies, 
are much more difficult to pinpoint. Only equili- 
brium reactions are necessary for steady-state, 
strictly diffusion-control flames. But to simulate 
the transient processes of ignition and flame 
propagation, kinetic reactions must be used, other- 
wise it was found that the flame propagation speeds 
will be unrealistically large. In its present con- 
figuration, the chemistry model in ARICC consists of 
13 kinetic reactions involving 8 chemical species, 
of which 9 are elementary reactions and 4 are third- 
body reactions. The elementary reactions are the 
same as those used by Westbrook (Ref. 9). These 
reactions are summarized in Table 4. The adequacy 
or redundancy of these reactions for various situa- 
tions are still being investigated. 

Boundary Conditions 


To complete this brief description of ARICC, 
mention should also be made of the boundary condi- 
tions currently in effect. For the purpose of 
modeling coaxial injector combustion chambers, ARICC 
is run in its axisymmetric mode with specified in- 
flow conditions. The exit plane, which in most 
cases is a sonic nozzle throat, is treated as a spe- 
cified pressure boundary. Other outflow conditions 
such as a zero-velocity gradient or enforced sonic 
conditions can also be specified. On the solid 
walls or internal solid boundaries, the boundary 
layer is approximated with a law-of-the-wall treat- 
ment. Sometimes a simple no-slip boundary condition 
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Table 4. List of Hydrogen/Oxygen Kinetic 
Reactions Currently Used in ARICC 


1. H + 0 2 

— * 

O + OH 

2. H 2 + O 



H + OH 1 

3. H 2 0 + 0 

k 

OH + OH 

4. H 2 0 + H 

, 

H 2 + OH 

5. H2O2 + OH 


h 2 o + ho 2 

6. H0 2 + H 


2 OH 

7. H0 2 + H 


H2 + 02 

8. H2O2 + O2 


2 H0 2 

9. H 2 0 2 + H 


ho 2 + h 2 

| THIRD BODY REACTIONS: j 

10. H + 0 2 + H 2 


ho 2 + h 2 

11. H + 20 2 

k 

H02 + 02 

12. H + 0 2 + H 2 0 

— - 

O 

CM 

X 

+ 

CM 

O 

X 

13. H2O2 + H2O2 


20H + H 2 0 2 1 


is preferred in the case of complicated surface con- 
tours . Details of the law-of-the-wall treatment are 
also discussed in Ref. 1. 


Examples 

Several cases of test runs are presented here to 
illustrate the qualitative features of the ARICC code. 

Propagation of Flame Front 

Figure 5 shows the model of a typical coaxial in- 
jector and combustion chamber. To check the one- 
dimensional flame propagation speed, the chamber is 
sealed off and filled with a premixed oxygen/hydrogen 
gas of 1:1 mixture ratio at 2000 psi and 298 K. At 
time 0, the entire upper row of cells is injected 
with a sufficient amount of energy to ignite the mix- 
ture, and the plane flame "sheet" propagates downward. 
It is interesting to note that a pressure wave 
travels at shock speed ahead of the flame front, 
which travels at the much lower speed of heat wave or 
concentration wave propagation. The flame front 
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FLAME < v tt ~ 120 lp,) 



t = 0.72S mwc 


Fig. 5. Propagation of Plain Flame Front 

speed is estimated to be approximately 120 ft/sec, com- 
pared with the theoretical estimate (through SSME data 
correlation) of 114 ft/sec. Also at later times, it 
is noticeable that the planeness of flame front starts 
to deteriorate, partly because of the presence of the 
injector at the axis and partly because of the compu- 
tational truncation errors toward the axis of symmetry. 

Transient Ignition in Prefilled Chamber 

In the second example, a similar combustion cham- 
ber has been prefilled with pure oxygen gas at 15 psi 
at 298 K. At time zero, a coaxial stream of gaseous 
hydrogen and oxygen at a mixture ratio of 1 starts 
coming in through the injector. Near the injector en- 
trance and slightly off the axis, a "spark" is ignited. 
The flow exits at the nozzle into atmospheric pressure. 
Figure 6 is a sequence of transient contour plots of 
water mass fractions, of hydrogen mass fractions, and 
of isotherms, and Fig. 7 is the velocity plot of the 
case early after startup showing the axisymmetric jet 
profile of the incoming hydrogen. The propagation of 
the major flame front, as diffused hydrogen reacts with 
the abundant oxygen, is evident. In this case, the 
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SECONDS 


A. H^O Mass Fraction 



0.0008 0.0012 0.0019 0.0030 0.0058 0.007® 

SECONDS 


B. OH Mass Fraction 



SECONDS 

C. Density 

Fig. 6. Transient Ignition in an Oxygen- 
Prefilled Chamber 
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VELOCITY VECTORS 
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Fig. 7. Velocity Vector Plot of Incoming 
Coaxial Jet Stream 

flame front is clearly defined by the sharp tempera- 
ture gradients and the presence of combusted products, 
such as water, with relatively uniform conditions in 
the combusted region behind. In addition to the major 
flame front, however, a smaller secondary and more 
permanent flame is being established right at injector 
entrance. This is the steady-state flame that is of 
more ultimate interest. At MR = 1, it is a fuel-rich 
flame, and both its size and intensity is dwarfed by 
the transient flame front propagating outward into the 
prefilled chamber. Yet as the transient flame propa- 
gates outward, its speed slows down, and it takes a 
long time (more than 20 msec) for it to consume all 
the prefill oxygen and burn Itself out. For this rea- 
son, another case is run in an attempt to simulate the 
steady state in a reasonably short time. 

Steady-State Hydrogen/Oxygen Flame 

In the third example, the same combustion chamber 
is now prefilled with the combustion products of the 
incoming hydrogen/oxygen stream at the estimated 
steady-state chamber pressure of 98 psi and tempera- 
ture of 797 K. The gaseous coaxial streams, this time 
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at the mixture ratio of 0.5, again comes in from the 
top. After approximately 3 msec, the "steady state" 
is well established. The isotherms, the H 0 mass frac- 
tion contours, and the OH mass fraction contours are 
shown in Fig. 8. The presence of OH is probably the 
indicator of where the chemically active region is. In 
this case, however, difficulty has been experienced in 
maintaining the flame. As soon as the spark is turned 
off, the flame has a tendency to quench out. Indeed, 
in comparable experiments, sometimes, though not con- 
sistently, the same problem of flame ignition has been 
experienced. Thus, it appears that the problem is re- 
lated to the details of the flame holding mechanism at 
the injector element tip, where the geometry is now 
only very coarsely modeled. A numerical model using 




Fig. 8. Steady-State Combustion of Gaseous 
Oxygen/Hydrogen Flame 

ARICC with very fine resolution around the injector 
should provide interesting insight into the exact 
nature of the flame holding mechanism there. 
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Liquid Jet Injection and Atomization 


In the final example, the unique capability of 
ARICC to handle two-phase flow and the creation of 
droplets from a liquid jet is demonstrated. Instead 
of gaseous oxygen, a stream of liquid oxygen (density 
= 69.97 lb/ft , temperature = 94.4 K) comes in through 
the center of the injector element. Gaseous hydrogen 
at atmospheric temperature comes in through the outer 
annulus at the high speed of 278 ft/sec (compared with 
4.4 ft/sec for the oxygen) to achieve a mixture ratio 
of one. The chamber is prefilled with hydrogen at 
100 psi and atmospheric temperature. Figure 9 is a 
sequence of plots showing the atomization process, and 

+ l 

» = 0.20 ms 0.30 ms 0.40 ms 

# PART. 183 276 348 


Fig. 9. Atomization Sequence of a 
Coaxial LOX Jet 

Fig. 10 is a closeup view of the evolving shape of the 
liquid jet. At this range of relative .velocities, 
droplets of the order of 10 microns are typically 
formed, mostly right in the recessed "cup" region 
shortly after injection. Note that the number of par- 
ticles indicated in the figure are numerical droplet 
groups, with each group representing a certain number 
of actual physical particles depending on the local 
stripping rate at time of creation according to the 



0.50 ms 
493 


0.53 ms 
518 
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t = 0.40 ms t = 0.50 ms 


Fig. 10. Closeup Details of an Evolving 
Liquid (LOX) Jet 

CICM atomization model. At the beginning of the run, 
the liquid jet column is given an arbitrary convenient 
shape (conforming to grid lines) with kinks on its 
surface. Toward the end of the run, at approximately 
0.4 msec, it can be seen that the fluid dynamic forces 
have elongated the liquid column and smoothed out its 
surface. Though no experimental data are available 
now for a direct comparison, the simulated jet size 
and shape and the spray cone angle of the droplets 
(approximately 25 degrees) are all intuitively correct. 

Conclusions and Recommendations 


The ARICC code has shown that with present-day 
computers and numerical methodology, it is now possi- 
ble to construct a complete model of the .combustion 
process inside a liquid rocket engine. The develop- 
ment process points to many areas for further work, 
not the least of which is the need for experimental 
correlations to be fed into the many physical models 
that still need to be described empirically. Much spe- 
cific measurements need to be done, using nonintrusive 
flow diagnostic techniques and spray measuring tech- 
niques, before the physical models can be confidently 
anchored and refined. 
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On the other hand, the availability of a tool like 
ARICC also opens up many potential applications. Al- 
ready mentioned is the use of ARICC as a comprehensive 
simulator for injector/combustor hardware and other com- 
bustion devices. Secondly, ARICC can be used as an en- 1 
gine performance analyzer, especially when small but 
powerful dedicated computers are available. In this ap- 
plication, ARICC combines the diverse capability of and 
will eventually supercede the earlier generation codes 
stated previously (Ref. 2-6). The third and most prom- 
ising area of application is ARICC as a research tool. 
While various physical models are now installed in 
ARICC, it can be turned around and used as a test bed 
for various upgraded models as new theories or data be- 
come available. Included in this category are areas 
for computational fluid dynamics research, such as tur- 
bulence models, algorithm efficiency improvements, and 
advanced graphics. Droplet dynamics research includes 
spray characteristics, droplet breakup/coagulation 
criteria, droplet deformation effects, atomization 
mechanism, and suspension flow or slurry flow. Thirdly, 
and most importantly, flame research using ARICC can be 
very fruitful in furthering understanding of diffusion 
flames, detonations, deflagrations, flame front propa- 
gation, ignition criteria, flame stability, and kinetic 
chemistry models. Finally, it should be remembered 
that ARICC is a two-dimensional/axisymmetric code. 

While already extremely valuable, most real-life hard- 
ware would require a fully three-dimensional model to 
study. Thus, the development of a three-dimensional 


counterpart to ARICC would be the logical next step to 


pursue. 
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